The hot non-perturbative gluon plasma is an almost ideal colored liquid 
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We study properties of a gluon plasma above the critical temperature T c in a generalized quasi- 
particle approach with a Lorentz spectral function. The model parameters are determined by a 
fit of the entropy s to lattice QCD data. The effective degrees of freedom are found to be rather 
heavy and of a sizeable width. With the spectral width being closely related to the interaction rate, 
we find a large effective cross section, which is comparable to the typical distance squared of the 
quasiparticles. This suggests that the system should be viewed as a liquid as also indicated by an 
estimate of the plasma parameter F. Furthermore, within the quasiparticle approach we find a very 
low viscosity to entropy ratio, r]/s ~ 0.2 for T > 1.05T C , supporting the recent conjecture of an 
almost ideal quark-gluon liquid seen at RHIC. 

PACS numbers: 12.38Mh, 25.75.-q 



The formation of a quark-gluon plasma (QGP) and its 
transition to interacting hadronic matter - as occurred 
in the early universe - has motivated a large community 
for more than two decades (cf. [1] and Refs. therein). 
However, the complexity of the dynamics in ultrarela- 
tivistic nucleus-nucleus collisions - producing high den- 
sity matter for short time scales - has not been fully un- 
ravelled and many properties of the new phase are still 
under debate [2]. In central Au+Au collisions energy 
densities are reached at the Relativistic Heavy Ion col- 
lider (RHIC) that are far above the critical energy den- 
sity e c ~ 1 GeV/fm 3 for a phase transition to a QGP as 
expected from lattice QCD calculations [3]. A strong ra- 
dial expansion and elliptic flow of hadrons, furthermore, 
point towards an early generation of pressure and a high 
interaction rate in the 'new phase' [4]. 

The latter observables are severely underestimated in 
conventional string/hadron transport models [5-7], how- 
ever hydrodynamical approaches do quite well in describ- 
ing (at midrapidity) the collective properties of the sys- 
tem for low and moderate transverse momenta [8]. The 
picture thus emerges that the medium created in ultra- 
relativistic nucleus- nucleus collisions for a couple of fm/c 
interacts more strongly than hadron/string matter, and 
it exhibits collective properties that resemble those of a 
liquid of low shear viscosity r\ [9]. In fact, viscous hy- 
drodynamical calculations indicate a very low viscosity 
to entropy ratio, rj/a ~ 0.1. ..0.2 [10]. This picture is 
substantially different from the expectation of a weakly 
coupled colored plasma. There is a variety of models that 
address the properties of this 'new matter'. It might be 
some kind of i) 'epoxy' [11], i.e. a system of resonant or 
bound gluonic states with large scattering length, ii) a 
system of chirally restored mesons, instanton molecules 
or equivalently giant collective modes [12], iii) a system 
of colored bound states of quarks q and gluons g, i.e. gq, 
qq, gg etc. [13]. 

In this letter we will provide quantitative arguments 
that strongly interacting matter in a certain tempera- 
ture range above T c is in a liquid phase. Our arguments 
are based on a generalized quasiparticle description of 



the system taking into account the spectral width 7 [14] 
in addition to the quasiparticle mass [15-19]. The model 
parameters are adjusted to non-perturbative results of 
lattice calculations. Since the width is closely related to 
the interaction rate we can then estimate relevant trans- 
port properties, such as the effective cross section and 
the shear viscosity, for temperatures near T c . Since lat- 
tice 'data' are more precise for quenched QCD we focus 
here on pure gluonic systems. We expect, however, simi- 
lar results for full QCD as argued below. The quantities 
we address are dominated by 'hard' momenta of the or- 
der of the temperature T. Accordingly, the quasiparticle 
properties we are interested in are related to the gluon 
propagator at hard momentum scales. In this line we 
take into account d g = 2(N 2 — 1) = 16 transverse gluons, 
and neglect Landau-damping contributions as well as the 
collective longitudinal modes whose spectral strength is 
suppressed for larger momenta [15, 16]. 

In order to adjust the quasiparticle properties we 
first consider thermodynamic bulk properties within the 
^-derivable formalism [20], which yields consistent re- 
summed approximations [21]. To leading-loop order, 
which is expedient for large coupling as argued in [14], 
the entropy follows directly from the quasiparticle prop- 
agator, cf. Ref. [15], 

(i) 

where n(uj) = (exp(u>/T) — 1) 1 denotes the Bose distri- 
bution function. We note that in the context of Fermi 
liquid theory a corresponding approximation is called the 
dynamical quasiparticle (dqp) entropy [22]. In principle, 
the resummed propagator A = (uj 2 — p 2 — n) _1 is to be 
calculated from a 1-loop Schwinger-Dyson equation. To 
proceed at this point, however, we use a physically mo- 
tivated Ansatz, assuming a Lorentzian spectral function, 

P{UJ) = E {{u-E) 2 +1 2 ~ {u J + E) 2 + 1 2 ) ■ (2) 
With the convention E 2 (p) = p 2 + M 2 — j 2 , the param- 
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eters M 2 and 7 are directly related to the real and imag- 
inary parts of the corresponding (retarded) self-energy, 
II = M 2 — 2ijoj. At this point we emphasize that the 
entropy functional (1) is not restricted to 'strict' quasi- 
particles, i.e., 7 need not to be small compared to the 
typical energy. Following the models [18] we parame- 
terize the quasiparticle mass in the gauge invariant and 
momentum independent asymptotic form 



M 



Nc 
6 



g 2 T 2 , 



with N c = 3 and the running coupling 



g 2 (T) 
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(3) 
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which permits an enhancement near T c [18, 23]. Likewise, 
we parameterize the width in the form 7 ~ g 2 T In g~ x [24] 
or, equivalently, in terms of M [14], 



7 4tt T 2 n (Al/T) 2 



(5) 



where c is related to the magnetic sector of QCD. We 
note that M and 7, as parameterizations of the complex- 
valued self-energy II at the relevant large momenta near 
the 'mass shell', are not related by a dispersion relation. 

In the upper part of Fig. 1 we compare the lattice 
results [25] to the quasiparticle entropy; the fitted pa- 
rameters are A = 2.42, T s = 0.46 and c = 14.4. We 
also display the interaction measure e — 3p (the pressure 
p and the energy density e are evaluated by thermody- 
namic relations, cf. [18]), which is particularly sensitive 
to interaction effects. We emphasize that the remarkable 
agreement with the lattice data is non-trivial because the 
functional relation between 7 and M is fixed, cf. (5). 

The adjusted quasiparticle mass and width are dis- 
played in the lower part of Fig. 1. The quasiparticles 
are rather heavy in line with direct lattice calculations 
[26]. For T ^5 1.05T C the width is sizeable, reaching more 
than 50% of the mass in a large temperature range. The 
picture of the strongly interacting plasma is, thus, a sys- 
tem of massive excitations with a large collisional width 
or short mean-free path - opposite to the original concept 
of narrow quasiparticles. Near T c , however, the width is 
close to zero. Although we have parameterized 7 in the 
'perturbative' form (5) we expect that the inferred tem- 
perature dependence is generic: near T c the width has to 
be small due to the small entropy at T m T c [14]. This 
is in line with a critical slowing down near a phase tran- 
sition. Away from T c one has to expect a large width 
due to the strong coupling and increasing reaction rates. 
The physical processes contributing to the width are then 
99 99 scatterings as well as splitting and fusion reac- 
tions gg <-> g or gg «-> ggg etc. Summing up the elastic 
and inelastic channel (and neglecting Bose-cnhancement 
for the final states) we end up with the total binary re- 
action rate 
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FIG. 1: The entropy s and the interaction measure w = e — 3p, 
in units of the Stefan-Boltzmann limits sq and po, from our 
quasiparticle model in comparison to lattice calculations [25]. 
The lower part shows the adjusted mass M and width 7. 



= (a) Tr Pl Tr P2 2^f\(s, P 2 , P 2 ) =: (a) 7 2 , (6) 

where A(x, y,z) = (x — y — z) 2 — Ayz and s = (Pi + P 2 ) 2 - 
In (6) we have introduced the shorthand notation 



Tr, 



duo d 3 p 
2^ (2tt) 3 



2up{oj)e(uj)Q(P 2 )n{uj) 



(7) 

for the thermally weighted trace over the quasiparticle 
degrees of freedom. The Q(P 2 ) function ensures that 
only time-like reaction processes are taken into account 
in Eq. (6). The interaction rate, on the other hand, is 
also related to the imaginary part of the self-energy; with 
a similar factorization as in Eq. (6), dN co u/dVdt — jN+. 
Here the particle density N + = Tr 1 is the time-like part 
of the integrated distribution function. The resulting ef- 
fective total cross section, 



(a) = 1 N + /I 2 



(8) 



is displayed in the top part of Fig. 2 as a function of T (us- 
ing T c — 0.26 GeV for quenched QCD). It rises from w 
at T = T c to about 20 mb at T w 1.1T" C , and drops again 
at higher temperatures. We note that similarly large val- 
ues for parton cross sections have been used in the phe- 
nomenological studies in Ref. [27]. These cross sections 
are larger by an order of magnitude than typical pertur- 
bative estimates for gg scattering (cf. the hatched band 
in the top part of Fig. 2) where the Debye mass mDebye is 
used as an infrared cutoff. However, for strongly coupled 
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plasmas the Debye mass may not be the proper regulator 
as pointed out by Thoma [28]. 
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FIG. 2: Upper part: The effective totaf cross section (<j), 
Eq. (8), in comparison to the perturbative estimate <t p qcd = 
97ra 2 /2mQ cbyc . 2nd part from top: The percolation measure 
k, Eq. (9); the critical value is n c — 1.18 [29]. 3rd part from 
top: The plasma parameter F, Eq. (10). Bottom part: The 
ratio of shear viscosity to entropy in our quasiparticle model 
in comparison to the lattice calculation [36] as well as to the 
next-to-leading log (nil) result [31]. Note that r\j s ~ 0.1 . . . 0.2 
was estimated from hydrodynamical fits to RHIC data [10]. 

A question of particular interest is the phase structure 
of the strongly interacting system, which in Refs. [9, 28] 
was surmised a liquid. In our quasiparticle approach we 
can address this question quantitatively by comparing 
the effective cross section to the typical distance squared 



of the quasiparticles to obtain information about critical 
clustering (percolation) [29]. If 



(a)N 



2/3 



(9) 



is lower than the critical percolation parameter k c ~ 1.18 
[29] the system is in the kinetic (dilute gas) regime while 
for k > k c percolation sets in and multi-particle interac- 
tions take over as characteristic for a liquid or solid with 
attractive interactions. As shown in the second part of 
Fig. 2, n is larger than k c for the temperature range 
1.05 . . .4T C . This suggests that the QCD plasma, up to 
rather large energy densities, is in a liquid/solid phase. 

In order to distinguish a liquid from, possibly, a solid 
phase we consider the plasma parameter [28] . Assuming 
equal (color) magnetic and electric energies we estimate 



T = 2 



N c g 2 



1 



47ri\T-V3 ( Tkm ) ' 



(10) 



with N being the full particle density (including the 
space-like contributions) , and the average kinetic energy 



(Thin) 



A+ 1 tr 



(11) 



Empirically it is well established for various attractive 
interactions that systems with 1 T iS 100 are in a liquid 
phase [30] . Since the 'critical' values may slightly vary for 
different systems we emphasize for the present case that 
the criteria r c = 1 for a gas-liquid phase transition and 
k = K c for percolation are met at the same temperature, 
cf. Fig. 2. Thus our results suggest that the gluon plasma 
is a liquid for temperatures between T c and 4T C . 

An important property of this gluon liquid is its shear 
viscosity rj. For weak coupling, it has been calculated in a 
transport approach from a Boltzmann equation to next- 
to-leading log order [31], if 111 ~ T 3 /(g 4 In g^ 1 ). However, 
for strongly coupled systems such an approach (assum- 
ing narrow quasiparticles) might be questionable and the 
Kubo formalism more appropriate instead [32] . Here the 
viscosity is evaluated from the slope of the Fourier trans- 
form of the spectral function ([Tij(x),Tij(y)]) for u> — > 0, 
where denotes the traceless part of the spatial stress 
tensor. At 1-loop order (and neglecting the longitudinal 
contributions) this corresponds directly to our quasipar- 
ticle picture which yields, cf. [33], 



dqp 



d g fdtu d 3 p dn 
2^ (2tt) 3 Ouj 



60 



P 2 {lo) [7uj 4 - lOu V + 7p 4 ] 



(12) 

Note that in the perturbative limit Eq. (12) (~ T 4 /-f ~ 
T 3 /(g 2 \ng- 1 )) does not approach rj nl1 - T 3 /\g 4 lug' 1 ). 
This comes about as follows: The shear viscosity rj is in- 
versely proportional to the transport cross section otrans 
in which the total scattering rate is weighted by (1— cos 9) 
with denoting the scattering angle [34]. In the weak 
coupling limit the gluon scattering is strongly forward 
peaked with a low gain in transverse momentum. This 
implies that many scatterings, corresponding to ladder 
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diagrams in the Kubo formalism, have to be resummcd 
in order to achieve a significant transverse momentum 
deflection and transport cross section. However, in the 
case of strong coupling the transport mean-frcc-path 
^trans ~ 1/ (&tran S N + ) is expected to become compa- 
rable to the total mean- free-path A ~ l/(aN+), as taken 
into account by Eq. (12). Thus, while a resummation of 
ladder diagrams is beyond the scope of the quasiparti- 
cle model, it should yield a useful approximation near T c 
[35]. 

In the bottom part of Fig. 2 we display the ratio of 
the shear viscosity to entropy. Although the lattice re- 
sults [36] still have large uncertainties, they are distinctly 
smaller than the (extrapolation of the) next-to-leading 
log result [31]. For T ^ 1.05T C our quasiparticle result is 
almost constant, ri dqp / s s=a 0.2 - in good agreement with 
the estimate r\j s « 0.1. ..0.2 from hydrodynamical fits to 
RHIC data [10]. The fast increase near T c , which is re- 
lated to the characteristic temperature dependence of 7 
and thus also of (a), can be seen as a precursor of the 
phase transition. We mention that the conjectured lower 
limit [37], rj/s > l/(4-7r), is approximately obtained for 
T « 1.1T C when calculating 77 in a non-relativistic limit of 
massive quasiparticlcs and assuming isotropic transport 
cross sections [35] . Since utrans < cr this estimate should 
give a lower bound for the viscosity. In conclusion, the 
quasiparticle model gives the picture of an almost ideal 
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The extension of the quasiparticle approach to the 
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the gg <-> g, gg <-> gg, gg <-> ggg and gg <-> gggg processes 
as also suggested in [39]. After approximate thermaliza- 
tion - of order 0.5 to 1 fm/c - the system behaves like an 
almost ideal massive colored parton liquid and exhibits a 
large pressure. This early pressure is responsible for the 
transverse flow of hadrons, and the large cross sections 
result in an almost complete suppression of far-side jets 
in central collisions. 
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